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Abstract 

In this paper, we consider massless Dirac fields propagating in the outer region of de Sitter- 
Reissner-Nordstrom black holes. We show that the metric of such black holes is uniquely determined 
by the partial knowledge of the corresponding scattering matrix S(X) at a fixed energy A / 0. More 
precisely, we consider the partial wave scattering matrices S(X, n) (here A 7^ is the fixed energy and 
n G N* denotes the angular momentum) defined as the restrictions of the full scattering matrix on 
a well chosen basis of spin-weighted spherical harmonics. We prove that the mass M, the square of 
the charge Q 2 and the cosmological constant A of a dS-RN black hole (and thus its metric) can be 
uniquely determined from the knowledge of either the transmission coefficients T(X,n), or the reflex- 
ion coefficients R(X, n) (resp. L(X, n)), for all n £ C where £ is a subset of N* that satisfies the Miintz 
condition ^2 nec — = +00. Our main tool consists in complexifying the angular momentum n and in 
studying the analytic properties of the "unphysical" scattering matrix S(X, z) in the complex variable 
z. We show in particular that the quantities T ^ ^ , and belong to the Nevanlinna class 

in the region {z £ C, Re(z) > 0} for which we have analytic uniqueness theorems at our disposal. 
Eventually, as a by-product of our method, we obtain reconstrution formulae for the surface gravities 
of the event and cosmological horizons of the black hole which have an important physical meaning 
in the Hawking effect. 
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1 Introduction 

Black hole spacetimes are among the most fascinating objects whose existence is predicted by Einstein's 
General Relativity theory and have attracted most attention in the last decades. From the theoretical 
point of view they are simple systems. The only needed parameters for a full description are the mass, the 
electric charge, the angular momentum and possibly the cosmological constant of the black hole raising 
the natural issue of determining them. From the astrophysical point of view however, black holes are 
objects eminently difficult to grasp since they are, by essence, invisible. Only by indirect means can we 
study some of their properties and... find actual evidence for their existence! A fruitful approach to 
better understand their properties consists in studying how black holes interact with their environment. 
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In particular, it is now well known established that much can be learned by observing how incoming waves 
are scattered off a black hole. We refer for instance to [21 [3j H21 US [H3 [22l 022 [33] where direct scattering 
theories for various waves have been obtained, to [U 03 [211 EH] for an application of the previous results 
to the study of the Hawking effect and to [6j [17] for an analysis of the superradiance phenomenon. In 
this paper, we follow this general strategy and address the problem of identifying the metric of a black 
hole by observing how incoming waves with a given energy A propagate and scatter at late times. This 
information is encoded in the scattering matrix S(X) introduced below. More specifically, we shall focus 
here on the special case of de Sitter- Reissner- Nordstrom black holes and we shall show that the parameters 
(and thus the metric) of such black holes can be uniquely recovered from the partial knowledge of the 
scattering matrix S(X) at a fixed energy A ^ 0. This is a continuation of our previous works [13l [H] in 
which similar questions were addressed and solved from inverse scattering experiments at high energies. 



1.1 de Sitter-Reissner-Nordstom black holes 

De Sitter-Reissner-Nordstrom (dS-RN) black holes are spherically symmetric electrically charged exact 
solutions of the Einstein-Maxwell equations. In Schwarschild coordinates, the exterior region of a dS-RN 
black hole is described by the four-dimensional manifold A4 = M. t x]r_, r + [ r xS 2 equipped with the 
lorentzian metric 

g = F{r)dt 2 - F{r)- X dr 2 - r 2 (d6 2 + sin 2 6 dtp 2 ) , (1.1) 

where 

2M Q 2 A 2 , . 

F(r) = l + *- r 2 . 1.2 

The constants M > 0, Q € K appearing in (|1.2p are interpreted as the mass and the electric charge of the 
black hole and A > is the cosmological constant of the universe. We assume here that the function F(r) 
has three simple positive roots < r c < r_ < r + and a negative one r n < 0. This is always achieved if 
we suppose for instance that Q 2 < §Af 2 and that AM 2 be small enough (see [2B])- The sphere {r = r c } 
is called the Cauchy horizon whereas the spheres {r = r_} and {r = r + } are the event and cosmological 
horizons respectively. We shall only consider the exterior region of the black hole, that is the region 
{r_ < r < r + } lying between the event and cosmological horizons. Note that the function F is positive 
there. 

The point of view implicitely adopted throughout this work is that of static observers located far 
from the event and cosmological horizons of the black hole. We think typically of a telescope on earth 
aiming at the black hole or at the cosmological horizon. We understand these observers as living on world 
lines {r = ro} with r_ << ro << r + . The variable t corresponds to their true perception of time. The 
event and cosmological horizons which appear as singularities of the metric (jl.ll) are in fact due to our 
particular choice of coordinates. Using appropriate coordinates system, these horizons can be understood 
as regular null hypersurfaces that can be crossed one way but would require speeds greater than that of 
light to be crossed the other way. From the point of view of our observers however, these horizons are 
thus the boundaries of the observable world. This can be more easily understood if we notice that the 
event and cosmological horizons are in fact never reached in a finite time t by incoming and outgoing 
radial null geodesies, the trajectories followed by classical light-rays aimed radially at the black hole or at 
the cosmological horizon. Both horizons are thus perceived as asymptotic regions by our static observers. 

Instead of working with the radial variable r, we make the choice to describe the exterior region of 
the black hole by using the Regge- Wheeler (RW) radial variable which is more natural when studying 
the scattering properties of any fields. The RW variable x is defined implicitely by *p = F~ l (r), or 
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explicitely by 



x = - — ln(r - r n ) + - — ln(r - r c ) + - — ln(r - r_) + - — ln(r + - r) + c, (1.3) 
where c is any constant of integration and the quantities Kj, j — n, c, — , + are defined by 

«n = ^'K), «c = ^(r c ), «_ = ^(r-), k + = \F'{r + ). (1.4) 

The constants k_ > and k + < are called the surface gravities of the event and cosmological horizons 
respectively. Note from (11.31) that the event and cosmological horizons {r = r±} are pushed away to the 
infinities {x — ±00} using the RW variable. Let us also emphasize that the incoming and outgoing null 
radial geodesies become straight lines {x = ±i} in this new coordinates system, a fact that provides a 
natural manner to define the scattering data simply by mimicking the usual definitions in Minkowski- 
spacetime. At last, note the presence of a constant of integration c in the definition of x. We shall 
comment on this constant and its consequences on our definition of the scattering matrix below. 



1.2 The scattering matrix and statement of the result 

As waves, we consider massless Dirac fields propagating in the exterior region of a dS-RN black hole. 
We refer to [301 [S3] for a detailed study of this equation in this background including a complete time- 
dependent scattering theory. We shall use the expression of the equation obtained in these papers as the 
starting point of our study. Thus the considered massless Dirac fields are represented by 2 components 
spinors ip belonging to the Hilbert space L 2 (R x S 2 ; C 2 ) which satisfy the evolution equation 

idtij = (r 1 D w + a(x)D S 2^, (1.5) 

where T 1 = diag(l, —1), D x = —id x and Dg2 denotes the Dirac operator on S 2 . Here, the potential a(x) 
takes the form 

¥) = - — 7—\ — ' ( L6 ) 

r(x) 

and thus contains all the information of the metric through the function F. In the variable x, it will 
be shown to have the following asymptotics a{x) ~ a±e K±x , x — > ±00 where a± are fixed constants 
depending on the parameters of the black hole. The equation (|1.5p is clearly spherically symmetric and 
in consequence can be separated. The stationary scattering is thus governed by a countable family of 
one-dimensional stationnary Dirac equations of the following form 

T 1 A, - (l + i)a(x)r 2 ] ip(x, A, = A tP{x, A, /), (1.7) 

restrictions of the full stationary equation to a well chosen basis of spin-weighted spherical harmonics 
(indexed here by I = |, |, ...) invariant for the full equation. Here T 1 and T 2 are usual 2x2 Dirac 
matrices satisfying the anticommutation relations T'T 3 + r j T l = A is the energy of the considered 
waves and (i + 5), I € i + N is called the angular momentum. For simplicity, we shall denote the angular 
momentum I + i by n. Hence the new parameter n runs over the integers N. 

As expected thanks to our choice of variable x, the stationary equation (jl.7j) is a classical one- 
dimensional massless Dirac equation in flat spacetime perturbed by an exponentially decreasing matrix- 
valued potential in which the angular momentum n plays the role of a coupling constant. Complete 



3 



stationary scattering theories have been obtained for this type of equation for instance in p] [20l [23] . 
Following the approach used in pQ, we can thus define in the usual way the scattering matrix S(X,n) 
in terms of stationary solutions with prescribed asymptotics at infinity, called Jost solutions. These are 
2x2 matrix-valued functions Fl and Fr solutions of (|1.7[) having the asymptotics 



F L (x,X,n) = e lTlXx (I 2 +o(l)), x -> +00, 
F R (x,X,n) = e trlXx (I 2 +o(l)), x -> -00. 



The Jost solutions will be shown to be fundamental matrices of (|1.7|) . There exists thus a 2 x 2 matrix 
Al(A, n) depending only on the energy A and the angular momentum n such that the Jost functions are 
connected by 

F L (x,X,n) = F R (x,X,n)A L (X,n). 

The coefficients of the matrix encode all the scattering information of equation (|1.7[) . In particular, 
using the notation 

a L1 (X,n) a L2 (X,n) 
a L3 (X,n) a L4 (X,n) 

the partial wave scattering matrix S(X,n) is then defined by 



A L (X,n) 



(1.8) 



S(X,n) 

where 



T(A, n) R{X,n) 
L(X,n) T(X, n) 



(1.9) 



tv\ ^ ^ \ a L2(X,n) a L3 (X,n) 

T(A,n) = a L1 (X,n), R{X,n) = — -, L{X,n) = — -. (1.10) 

oli(A, n) a,Li{X,n) 

The quantities T and R, L are called the transmission and reflection coefficients respectively. The former 
measures the part of a signal transmitted from an horizon to the other in a scattering process whereas 
the latters measure the part of a signal reflected from an horizon to itself (event horizon for L and 
cosmological horizon for R j]]. At last, the scattering matrix S(X, n) will be shown to be a 2 x 2 unitary 
matrix. 

Roughly speaking the main result of this paper states that either the knowledge of the transmission 
coefficient T(A, n) or the knowledge of the reflection coefficients L(A, n) or i?(A, n) at a fixed energy A 7^ 
and "for almost all" n € N determines uniquely the mass M and the square of the charge Q 2 of the 
black hole as well as the cosmological constant A of the universe. More precisely, it suffices to know the 
transmission or reflection coefficients at a fixed energy A ^ on a subset £cN* that satisfies the Miintz 
condition Xme.c h = 00 m or der to prove the uniqueness of the parameters M, Q 2 , A. Since the data of 
the partial wave scattering matrices 5(A, n) for all n G N is equivalent to know the full scattering matrix 
5(A), we can rephrase our main result by: the partial knowledge of the scattering matrix S(X) at a fixed 
energy A 7^ determines uniquely the metric of a dS-RN black hole. 

Before entering in the description of the method used to proved the above uniqueness result, let us 
comment on its dependence with our choice of coordinates system. As already mentioned, the variable x 
is defined by (|1.3j) up to a constant of integration c. Our definition (|1.9|) - ()1.10|) of the scattering matrix 
turns out not to be invariant when we change the constant c in the definition of x. More precisely, if 
S(X, n) denotes the scattering matrix obtained for a given RW variable x, we can show that the scattering 
matrix S(X) obtained using the translated RW variable x = x + c is given by 

S{X, n) = e lFlXc S(X, n) e - lFlXc , (1.11) 



1 Whence the notations L for left reflection coefficient since the event horizon is located "on the left'" at x = —00 and R 
for right reflection coefficient since the cosmological horizon is located "on the right" at x = +00. 
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or written in components by 



" T(A, n) R(X,n) 
L(X, n) T(X, n) 

Since there is no natural - and better - way to fix the choice of the constant c in (| 1 .3[) , we must include the 
possibility to describe a dS-RN black hole by two different RW variables in the statement of our result. 
One way to make our result coordinate invariant is to identify the partial wave scattering matrices at 
a fixed energy A corresponding to all the possible choice of RW variables in the description of a given 
dS-RN black hole. In other words, we shall say that S(X, n) and S(X, n) are equal when {TjTT} or (|I.I2[) 
hold. 

Having this in mind, we state now the main uniqueness result of this paper. 

Theorem 1.1. Let (M, Q, A) and (M, Q, A) be the parameters of two dS-RN black holes. We denote by 
a{x) and a(x) the two corresponding potentials appearing in the Dirac equation il.5\) . We also denote by 
S(X, n) and S(X, n) the corresponding partial wave scattering matrices at a fixed energy A 7^ defined by 
il.9\) and hi. 10)) . Consider a subset C ofW that satisfies the Miintz condition X)ne£ \ = 00 an< ^ assume 
that there exists a constant c G K such that one of the following conditions holds: 

(i) T(A, n) = T(A, n), Vn G C, 
(ii) L{X,n)=er 2lXc L{X,n) 1 Vn G C, 
{Hi) R(X, n) = e MXc R(X, n), Vn G C. 

Then the potentials a and a coincide up to translation, i.e. there exists a constant a G M such that 

a(x) = a(x + a), G K. 

As a consequence we get 

M = M, Q 2 = Q 2 , A = A. 
Let us make several comments on this result. 

1) We emphasize that the uniqueness results in Thm 11.11 are in fact twofold. First, we prove that 
a positive exponentially decreasing potential a(x) satisfying (I2.5l) - (|2.6p for the Dirac equation (jl.5[) is 
uniquely determined (up to translation) from one of the assumptions (i)-(iii) of Thm fTTT] Then and only 
then do we use the particular expression (|1.6|) of this potential to show that the parameters of the black 
hole (and thus the metric) are uniquely determined. 

2) As a particular case of Thm 11.11 we see that the potential a{x) and thus the parameters of the 
black hole are uniquely determined (up to translation for the potential) by the full scattering matrix 
S(X) at a fixed energy A ^ 0. Theorem 11.11 is sharp in the sense that the full scattering matrix S(0) 
at the energy A = does not determine uniquely the potential and the parameters (see Remark l3.7l below) . 

3) In the case of nonzero energies, it is also natural to ask whether the Miintz condition Y) ne £ — = 00 
is necessary? On one hand, as regards the problem of uniquely determining the parameters of the black 
hole, it is likely that we could weaken this condition since the metric only depends on "three" parame- 
ters. On the other hand, as regards the problem of uniquely determining the potential a(x), the Miintz 
condition could be sharp. Indeed a similar inverse scattering problem for 3D Schrodinger operators 
with radial potentials has been already studied by A. G. Ramm in [33] and M. Horvatz [23]; Ramm 



T(X,n) 
3~ 2lAc Z(A,n) 



R(X,n) 



T(X,n) 



(1.12) 
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showed that the knowledge of a subset of the phase shifts Si, with ^2 leC r = oo, determines uniquely the 
potential; shortly after, Horvath proved the necessity of the Miintz condition in some classes of potentials. 

4) At last, let us say a few words on our uniqueness results from a more geometrical point of view. 
Notice first that the Dirac equation (| 1 . 5[) in the exterior region of a dS-RN black hole takes the same 
form as a Dirac equation on the manifold £ = R x x S 2 equipped with the riemanniann metric 

g a = dx 2 + a - 2 (x)(d8 2 + sin 2 8dtp 2 ), (1.13) 

where a(x) is any smooth positive function. If we assume moreover that the function a(x) has the 
asymptotics (|2.5[) - (I2.6[) as it is the case in our model, then the riemanniann manifold (£, go) can be viewed 
as a spherically symmetric manifold having two ends {x = ±00} that are asymptotically hyperbolic. Hence 
our model fits the more general framework of asymptotically hyperbolic manifolds (AHM). In this setting, 
Thm II .ll states that metrics like (|1.13|) are uniquely determined (up to translations in x) from the partial 
knowledge of the scattering matrix S(X) - corresponding to Dirac waves - at a fixed energy A ^ 0. For 
more general AHM with no particular symmetry, some direct and inverse scattering results - for scalar 
waves - have been obtained by Joshi, Sa Barreto in [27] and by Sa Barreto in [36] (see also [26] and [8]). 
In |27] for instance, it is shown that the asymptotics of the metric of an AHM are uniquely determined 
(up to diffeomorphisms) by the scattering matrix S(X) at a fixed energy A off a countable subset of R. 
In [36], it is proved that the metric of an AHM is uniquely determined (up to diffeomorphisms) by the 
scattering matrix 5(A) for every A € R \ 0. 



1.3 Overview of the proof 

The main idea of this paper is to complexify the angular momentum n = I + i and study the analytic 
properties of the "unphysical" scattering coefficients T(A, z), L(X, z) and R(X, z) (or equivalently the 
functions aLj(X,z)) with respect to the variable z € C. The general idea to consider complex angular 
momentum originates in a paper by Regge |35] as a tool in the analysis of the scattering matrix of 
Schrodinger operators in R 3 with spherically symmetric potentials. We refer to [32] . chapter 13, and [9] 
for a detailed account of this approach. Applications to the study of inverse scattering problems for the 
same equation can be found in [TO] [TT] [34] . These last papers were the starting point of our work. 

The first step in our proof of Theorem 11.11 relies on uniqueness theorems for analytic functions. Let 
us define the Nevanlinna class A^(n + ) as the set of all analytic functions f(z) on the right half plane 
n+ = {zeC: Re{z) > 0} that satisfy 



sup / In 

0<r<l , 



1 



1 + re l f 



dip < GO, 



_!_ ! In x In x ^* 

where In (x) = < ^ ' 1 ' < q' Among other properties, it turns out that such functions are uniquely 

determined by their values on any subset £cN* that satisfies the Miintz condition X)ne£ n = 00 ( see 
[33] and [35]: chapter 15, for a more general statement). We shall use the Nevanlinna class and this 
uniqueness result as follows. 

Using explicit representations for the Jost functions as well as the unitarity of the scattering matrix, 
we first show that the coefficients aLj(\, z), j = 1, ..,4 in (11.8[) are entire functions of exponential type in 
the variable z that satisfy the bound: 

\a Lj (X,z)\<e A ^^, 



where A is the constant given by J R a(x)dx. From these estimates, we deduce easily that the functions 
ciLj(A, z) restricted to the right half plane II + belong to the Nevanlinna class 7V(II + ). Hence the preceding 
discussion allows us to conclude that the functions a,Lj(A,z) are completely determined by their values 
on any subset £cN* such that J2 n ec n = 00 ■ 

Since the true scattering data are the transmission coefficient T or the reflection coefficients L, R and 
not exactly the aLj, we need to work a bit more to get a useful uniqueness statement. Using mainly 
Hadamard's factorization theorem and the previous result, we show that in fact, the whole matrix Al (A, z) 
is uniquely determined from the values of one of the scattering coefficients T(A, n), L(X, n) or R(\, n) on 
any subset £cN* such that J2 n ec n ~ 00 on ly- 

The second step in our proof relies on precise asymptotics for the coefficients Oij(A, z) when the 
parameter z — > 00 for real values of z. To obtain these asymptotics, it is convenient to introduce a new 
radial variable X that has also the great interest to enlight the underlying structure of equation (|1.7[) . 
Following (101 ITT) , we define the variable X by the Liouville transformation 

X = I a(s)ds. (1.14) 

J — 00 

Note that X is well defined thanks to the exponential decay of a(x) at both horizons and runs over the 
interval (0,A) with A = J R a(s)ds. Let us denote by h(X) the inverse transformation of (| 1 . 14[) . We 
shall also use the notations Fl(X) and Fr(X) as a shorthand for the Jost functions Fi,(h(X), A, z) and 
Fn(h(X), A, z). The reason why we introduce such a variable lies in the observation that the components 
fLj{X) and fnj(X) of the Jost functions 



Fl{X) = 



hi(X) f L2 {X) 
f L3 (X) f Li (X) 



Fr(X) = 



fm(X) f m (X) 
fm{X) f m (X) 



satisfy then second order differential equations of the form 

f"(X) + q(X)f(X) = z 2 f(X), l£(0,4 (1.15) 
Here the potential q will be shown to have quadratic singularities at the boundaries and A, i.e. 
q(X)^^ 1 X^0, q(X)~ X^A, 

where uj± are two constants. We emphasize that the angular momentum or coupling constant z has now 
become the spectral parameter of this new equation. We shall show by a perturbative argument that 
the Jost functions fi,j(X) and fnj(X) can be written as small perturbations of certain modified Bessel 
functions^ From the well known asymptotics of the modified Bessel functions for large z, we then obtain 
precise asymptotics for the Jost functions when z — ¥ 00, z real, which in turn immediately yield the 
asymptotics of the coefficients Oij(A, z). Let us remark here that we could also obtain asymptotics of 
the scattering data for large values of z in the complex plane but we don't need such asymptotics in our 
proof. At last, we mention that singular Sturm-Liouville operators like (|1.15p have been studied in details 
by Freiling and Yurko in |18) including the solution of various inverse spectral problems. 

Let us now explain briefly how we prove our main theorem. Consider two dS-RN black holes with 
parameters M, Q, A and M,Q,A respectively. We shall use the notation Z and Z for all the relevant 



2 Note that the modified Bessel functions are solutions of the free equations obtained from 11.151 by replacing the potential 
q(X) by its asymptotics at the boundaries, i.e. f"{X) + ^rf(X) = z 2 f(X) when X -> and f"(X) + "+ . 2 f(X) = 

z 2 f(X) when X -> A. 
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scattering quantities relative to these black holes. Assume that one of the conditions in Thm 11.11 holds. 
From the previous uniqueness result, we conclude first that A L (X, z) = e~ 4r Ac i L (A, z)e* r Ac for all z e C. 
Second, we use a standard procedure in one-dimensional inverse spectral problem (see [TH]) together with 
the precise asymptotics for the functions Al(\,z) and Al(X,z) obtained previously to prove that there 
exists k € Z such that 

x + c + — J , Vx G K. 

Hence the potentials a and a coincide up to translations. This proves the first part of Thm fTTTl Using the 
particular expression of the functions a and a in terms of the parameters of the black holes, we deduce 
then that 



M = M , Q 2 = Q Z , A = A. 

Finally, as a by-product of our method, we also obtain explicit reconstruction formulae for the surface 
gravities k± from the reflection coefficients L(X,n) or R(X, n) when n — > oo. 

This paper is organised as follows. In Section [21 we recall all the direct scattering results of [T| 130] 133] 
useful for the later analysis. In Section [3] we put together all the analytical results concerning the 
scattering data At (A, z) and S(X, z). In particular, we show there the uniqueness theorems stated above. 
In Section^ we introduce the Liouville variable X and calculate the precise asymptotics of the scattering 
data Al(X,z) when z — > oo. In Section [5] we solve the inverse problem and prove our main Theorem 
ThmO 



2 The direct problem 

In this section, we first recall the expression of the Dirac equation in dS-RN black holes as well as the 
direct scattering theory obtained in (30] [33] . We then give an explicit stationary representation of the 
related scattering matrix following the approach and the notations used in pQ. 

As explained in the introduction, we describe the exterior region of a dS-RN black hole using the 
Regge- Wheeler variable x defined in (|1.3p . We thus work on the manifold B = K. t x E with E = x S 2 , 
equipped with the metric 

g = F(r){dt 2 -dx 2 )~r 2 duj 2 , 

where F is given by (II. 2|) and duj 2 = d9 2 +sin 2 9 dip 2 is the euclidean metric induced on S 2 . The manifold 
B is globally hyperbolic meaning that the foliation E t = {t} x E by the level hypersurfaces of the function 
t, is a foliation of B by Cauchy hypersurfaces (see [33] for a definition of global hyperbolicity and Cauchy 
hypersurfaces). In consequence, we can view the propagation of massless Dirac fields as an evolution 
equation in t on the spacelike hypersurface S, that is a cylindrical manifold having two distinct ends: 
{x = —oo} corresponding to the event horizon of the black hole and {x = +00} corresponding to the 
cosmological horizon. Let us recall and emphasize here the nature of the geometry - of asymptotically 
hyperbolic type - of the hypersurface E near the horizons. This peculiar geometry will be seen in the 
asymptotic behaviour of the potential a(x) below. 

2.1 Dirac equation and direct scattering results 

We use the form of the massless Dirac equation obtained in [30] [33] . The fields are represented by a 2- 
components spinor belonging to the Hilbert space H — L 2 (M.x S 2 ,dxdu>; C 2 ) and the evolution equation 
can be written under Hamiltonian form as 

id t ^ = H^, (2.1) 
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where the Hamiltonian H is given by 

H = T 1 D x + a(x)D S 2. (2.2) 



The symbol D x stands for — id x whereas D52 denotes the Dirac operator on S 2 which, in spherical 
coordinates, takes the form 

cot ( 

~2~ ' sin* 

The potential a is the scalar smooth function given in term of the metric (|l.l[) - (|1.2p by 



D S , = ^Y\d e + -^r 3 a v . (2.3) 



where r(x) is the inverse diffeomorphism of (|1.3|) . Finally, the matrices r x ,r 2 ,r 3 appearing in (|2.2[) and 
are usual 2x2 Dirac matrices that satisfy the anticommutation relations 



r i r j + r j r l = 2%, v», j = 1, 2, 3. 

We shall work with the following representations of the Dirac matrices 

r 1 = 





, r 2 = ( 


f 


i 




i \ 












-i ) 



We use now the spherical symmetry of the equation to simplify further the expression of the Hamil- 
tonian H. We introduce a basis of spin weighted spherical harmonics invariant through the action of 
the Hamiltonian that " diagonalize" the Dirac operator Dg2. We refer to I.M. Gel'Fand and Z.Y. Sapiro 
[19] for a detailed presentation of these generalized spherical harmonics and to [30, 33j for an application 
to our model. For each spinorial weight s, 2s E Z, we consider the family of spin-weighted spherical 
harmonics {Yj m = e m(p u l sm: I — \s\ G N, I — \m\ G N} which forms a Hilbert basis of L 2 (S 2 ,dui) and 
where the u l sm satisfy the following relations 



dul m m — s cos 9 , ., , ,1 , 

—7?r — a < m = -f + S i-S + lW-l,m> 

at) sm 

du[ m m — scos9 , 1 , 

-ST + — a — «L = -*[(i + a + l)(i-*)] 3 «i+i, m - 

o(t7 sm [/ 



We define the familly = m? Y| m ) where the indexes l,m belong to X = {(i,m), 2 — | G 

N, Z — \m\ G N}. Moreover, we define 02 as the following operation between two vectors of C 2 

Vu = (ltl,M 2 ),W = {vi,V 2 ), U® 2 V = {uiVl,U 2 V 2 )- 

Then the Hilbert space % can then be decomposed into the infinite direct sum 

H= Hi m , Hi m = L 2 (R x ;C 2 )® 2 Ft n . 

We shall henceforth identify Hi m and L 2 (R;C 2 ) as well as ipi m ® 2 F m an d i>i m - What's more, it is 
easy to check that the Hi m are let invariant through the action of H . Hence we obtain the orthogonal 
decomposition for the Hamiltonian H 

H= H lm , 

(Lm)eX 
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with 

H lm := H ]nim = T 1 ^ - (/ + \)a{x)T 2 . 

In what follows, it will be thus enough to consider the restriction of the Dirac equation (|2.ip to each 
Hilbert space T-Li m separatly and define there the partial-wave scattering matrices. The full scattering 
matrix will be then written as the orthogonal sum of these partial-wave scattering matrices. 

For later use, we precise here the asymptotics of the potential a in the RW variable defined by (|1.3|) . 

Lemma 2.1. 

a(x) = a ± e K±x + 0(e 3K±x ), x -> ±00, (2.5) 
a'{x) = a±K±e K±x + 0(e 3K±x ), x -> ±00, (2.6) 

where the constants k± are given by [1.4]) and a± by 



with 



C± = In 



Vr2^ e _ K±c±j (27) 



{r±-r n )z«n(r ± -r c )^(r + -r-) 2 ^ +c, (2.8) 



and c is the constant of integration from [1.3\) . 

Proof. Since the function F is naturally expressed using the r radial variable, we need to find first 
equivalents between the RW variable x and r when r — > r± or equivalently when x — ¥ ±00. From (jl.3[) . 
we have 

x = - — In \r — r±\ + C± + 0(\r — r±\), r— > r±, 
where C± are the constants given by (I2.8p . Hence we get 

e 2K± ^- c ^ =\r-r ± \(l + 0(\r-r ± \fj, r -> r±, (2.9) 
Taking the square root of (|2.9p . we also have 

e «±(x-C±) = ^|r-r±|(l + 0(|r-r±|)), r -> r±. (2.10) 
Now a simple calculation shows that 



a(x) = ^EE>. = ^ T ' 2K± ^\r-r±\ + 0(\r-r±\%), r -> r±. (2.11) 
r r± 

Hence f|2.5|) follows from (|2.10p and (|2.1ip . Similarly, we have 



a ' {x) = F{ r)±(yS\ =K± ^^^^ + 0(|r-r ± |t), r r ±J 
ar \ r I r± 

which together with ([2~T0|) yields (|2l)J) . □ 
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In particular, the potential a(x) is thus exponentially decreasing at both horizons. In consequence, the 
Hamiltonian H lm can be viewed as a very short-range perturbation of the free Hamiltonian Hq = T 1 D X 
on each spin weighted spherical harmonic. Using this fact, the spectral and scattering properties of the 
Hamiltonians H, H lm are established in an elementary way. The main scattering results obtained in 
[301 133) are summarized in the next proposition 

Proposition 2.2. (i) The Hamiltonians H lm and H are selfadjoint on Him and H with domains 
D(H lm ) = H^C 2 ) andD(H) = {ijj eH,f = ^im, frm € D{U lm ), £ {\\H lm i> lm \\ 2 + ||^ m || 2 ) < 

oo}. 

(ii) The Hamiltonians H lm ,H have no pure point and singular continuous spectra, i.e. 

a pp {H lm ),a pp (H) = 0, <J smg (H lm ),a smg (H) = 0. 
In other words, the spectra of H lm and H are purely absolutely continuous. 

(Hi) Denote by Hq = the selfadjoint operator acting on Him as well as on H. Then the wave 

operators Wj^ and W ± defined by 

W,t=s- Km e ltH ' m e~ ltH ° , W ± = s - lim e ltH e - ltH °, 

exist and are asymptotically complete. 

As a direct consequence of Proposition 12.21 we can define the partial wave scattering operators Si m 
by the usual formulae 

si m = (W+ ywf m , 

which are well-defined unitary operators on Hi m - Let us now introduce the unitary transform on Him 

(F i>i m )(\) = / e- iVlx *i> lm {x)dx. (2.12) 



2tt 

The transform F clearly diagonalizes the free Hamiltonian H — T l D x on Hi m . Note in passing that F 
also acts in a trivial way on H. Morover it is unitary and still diagonalizes H Q on H. By definition, the 
partial wave scattering matrices Si m (X) are simply the partial wave scattering operators Si m written in 
the energy representation of H given by (I2.12[) , that is for all 4n m £ £ 2 (Ra,C 2 ) we have 

SlmW&mW = (F Sl m FQ(j)lm)W- ( 2 -13) 

The partial wave scattering matrices Si m (X) are unitary 2x2 matrices for all A £ R. 

Remark 2.3. Using the spherical symmetry of the equation and the definitions of the wave operators in 
Proposition \2.2\ the full scattering operator S and the full scattering matrix 5(A) are unitary operators on 
H. and L 2 (S 2 , C 2 ) defined as the orthogonal sum of the partial wave scattering operators Si m and matrices 
Sim(X) respectively, i.e. 

s = (w + yw~= s lm , 

(Lm)eC 

S{\)<K\) = (FoSF$<f>)(\)= Si m {X)<f>i m {X), y<f> £ L 2 (M.\ x S 2 ,C 2 ). 

(l,m)eC 

Remark 2.4. Observe that the Hamiltonians H lm only depend on the angular momentum l + I £ 
From the above definitions, we see immediately that the wave operators Wi m , the scattering operators Si m 
and the scattering matrices Si m (X) only depend on the angular momentum I + h too. For simplicity we 
shall therefore denote by S(X,n), n = I + | £ N. the family of scattering matrices Si m (X). 
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2.2 A stationary representation of the scattering matrix 

In this section, we follow the approach of [T] and obtain an explicit stationary representation of the 
scattering matrices S(X, n) for a fixed energy A £ R and all angular momentum neN. Let us emphasize 
here that the results of this section still hold if we assume a £ L 1 (R) only. 

Let us consider first the stationary solutions of equation (|2.1|) restricted to each spin weighted spherical 
harmonic, i.e. the solutions of 

[r 1 ^ - na(x)T 2 }^ = Xtp, Vn € N. (2.14) 

Here we can think of ip in (|2.14[) as either a column vector of 2 entries, or as a 2 x 2 matrix. For A G M, 
we define the Jost solution from the left Fl(x, A,n) and the Jost solution from the right F R (x, X,n) as 
the 2x2 matrix solutions of (|2.14[) satisfying the following asymptotics 

F L (x,X,n) = e irlXx (I 2 + o(l)), x^ +oo, (2.15) 
F R (x, A, n) = e^^ih + oil)), x^ -oo. (2.16) 

From (|2.14l) . (|2 . 15[) and (|2.16|) . it is easy to see that such solutions (if there exist) must satisfy the integral 
equations 

r+oo 

F L (x,X,n)^e^ Xx -tnT 1 / e~ iTl x ^ y ~ x "> a(y)T 2 F L (y , A, n)dy, (2.17) 

J X 

F R (x, A, n) = e lTlXx + mV 1 f e - trlx ^a(y)T 2 F R (y, A, n)dy. (2.18) 



Since the potential a belongs to £ 1 (M), it follows that the integral equations (|2.17l) and (|2.18l) are uniquely 
solvable by iteration and that 

\\F L {x,X,n)\\ <e n L + °° a ^ ds , \\F R {x, A, n)\\ < e n ^ a(s)ds . 
Moreover we can prove 

Lemma 2.5. For A £ 1 and n € N, either of the Jost solutions Fl(x, A, n) and F R (x, X,n) forms a 
fundamental matrix of and has determinant equal to 1. Moreover, the following equalities hold 

F L (x, A, n)* T 1 F L (x, A, n) - T 1 , (2.19) 
F R (x,\n)*T l F L (x,X,n) = T\ (2.20) 

where * denotes the matrix conjugate transpose. 

Proof. See [I], Proposition 2.2. □ 

Since the Jost solutions are fundamental matrices of (|2.14l) . there exists a 2 x 2 matrix Al(X, n) 
such that Fl(x, A, n) — F R (x, A, n) Al(X, n). From (|2.16[) and (|2.17j) . we get the following expression for 
A L {X,n) 

A L (X, n)=I 2 - mT 1 f e- irlx y a (y)T 2 F L (y, A, n)dy. (2.21) 
Moreover, the matrix A^(A, n) satisfies the following equality (see [I], Proposition 2.2) 

A^(A,n)rM L (A,7i) = T 1 , VA € K, n £ N. (2.22) 
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Using the notation (jl.8|) . the equality (|2.22|) can be written in components as 



|a L i(A,n)| 2 - |a L3 (A,n)| 2 
|a L4 (A,n)| 2 - |a L2 (A,n)| 2 



= 1, 
_ = 1, 

«li(A, n)a i2 (A,n) - a L3 (X,n)a Li (X,n) = 0. 



(2.23) 



As mentioned in the introduction, the matrices A^(X, n) encode all the scattering information of equation 
(|2.14j) . In particular, it is shown in [I] that the scattering matrix S(X,n) defined in (|2.13l) has the 
representation 

" T(X, n) R(X,n) 
L(X,n) T(X,n) 



S(X,n) = 



where 



T(A, n) = a L1 (A, n), R(X, n) = — -, L(X, n) = — -. 

a L1 (X,n) a L1 {X,n) 



(2.24) 
(2.25) 



The unitarity of the scattering matrix S(X, n) leads to the following relations 
Lemma 2.6. For each A G K and neN, we have 



|T(A,n)| 2 + |i?(A,n)| 2 
|T(A,n)| 2 + |L(A,n)| 2 



(2.26) 



T{X,n)R{X 1 n) + L{X,n)T{X,n) = 0. 

Note that the relations (|2.26j) are also direct consequences of the relations (|2.23|) and definitions (|2.25l) . 

Remark 2.7. We finish this section analysing the influence of the constant of integration c used in the 
definition hi. Sty of the RW variable on the expression of the scattering matrix S(X,n). Assume thus that 
we describe the same dS-RN black hole using two RW variables x and x = x + c. We shall denote by Z 
and Z all the relevant quantities expressed using the variables x and x respectively. We follow the same 
procedure as above to define the scattering matrix S(X,n). Our goal is to find a relation between S(X,n) 
and S(X,n). We start from the stationary equation \2.11$ obtained using x. Hence we get 

[r 1 ^ - na{x)T 2 ]^ = Xtjj, 

where the potential a(x) is simply the translated by c of the potential a, i.e. a(x) = a(x — c). Therefore, 
the uniqueness of the Jost functions satisfying the asymptotics \2. 1 5\) and \2. 1 6\) yields 



F L (x,X,n) = F L (x - c,X,n)e 



Hence, it follows from (2.21]) that 



A L {X,n) = e~ lTlXc A L {X,n)e lTlXc 



Eventually, using ([ST. 



S(X,n) 



e- lTlXc S(X,n)e 



we conclude that 
r x Ac _ 



T(A, n) 
* 2lXc L(X,n) 



- 2lXc R(X,n) 
T(X,n) 



that is S7U\) . 
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3 Complexification of the angular momentum 



In this section, we allow the angular momentum to be complex. After studying the analytic properties 
of the Jost functions Fl(x, A, z), F R (x, A, z) and of the matrix Al(X, z) in the variable z £ C, we prove 
the uniqueness results mentioned in the introduction. 



3.1 Analytic properties of the Jost functions and matrix Al(X,z) 

Let us start with the Jost functions Fl(x, A, z) and F R (x, A, z). They are solutions of the stationary 
equation 

[r 1 ^ - za{x)T 2 ]i) = Xip, Vz G C. (3.1 



with the asymptotics (|2.15p and (|2.16p . By commodity, we introduce the Faddeev matrices Ml(x, A, z 
and M R (x, A, z) defined by 



M L (x,X,z) =F L (x,X,z)e- lFlXx , M R (x, A, z) = F R {x, A, z)e~ iFlXx , 



which thus satisfy the boundary conditions 



M L (x,X,z) = I 2 + o(l), x^+oo, 
M R (x,X,z) = I 2 + o(l), x — > — oo. 



We shall also use the notations in components 
M L (x,X,z) 



m L i(,x,X,z) m L2 (x,X,z) 
m L3 (x,X,z) m Li (x,X,z) 



M R (x,X,z) 



m m (x,X,z) m R2 (x,X,z) 
m R3 (x,X,z) m RA (x,X,z) 



From (|2.17|) and (|2.18|) . the Faddeev matrices satisfy the integral equations 



M L (x, A, n)=h- izT 1 / e- irl ^-^a(y)T 2 M L (y, A, z)e irlx ^dy, 



M R (x, A, n) =h + izT 1 / e-^ lx ^a{y)T 2 M R {y, A, n)e lTlx ^dy. 



Iterating ()3.6|) and (|3.7p once, we get the uncoupled systems 

p + OO Z' + OO 



m L i(x,X,z) 
m L2 (x, A, z) 
m L3 (x, A, z) 
rriL4{x, A, z) 



= l + z< 



^~y) a (y)a(t)m L1 (t,X,z)dtdy, 



OO p-\-OQ 



(3.2 



(3.3 
(3.4 



(3.5 

(3.6 
(3.7 

(3- 



-iz / e- 2lX{y - x) a(y)dy + z 2 / / e - 2lX(y -^a(y)a(t)m L2 (t, A, z)dtdj$.9 

J x J x J y 

iz e 2lX{y - x) a(y)dy + z 2 / e 2tX( - y - x ^a(y)a(t)m L3 (t, A, z)dtdy, (3.10 

J x J x J y 

f + OO /* + oo 



l + z z 



- 2tX ^- y ^a{y)a{t)m L4 (t, A, z)dtdy, 



(3.11 
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and 



/ e- 2lX{y - t) a{y)a(t)m R1 (t,X,z)dtdy, (3.12) 

-oo J — oo 

m R2 (x,X,z) = iz / e 2tXi - x ^a{y)dy + z 2 / / e 2lX( - x ~ y) a(y)a(t)m R2 (t, A, z)dtdy, (3.13) 



oo •/ — oo 



m fl ,3(x,A,z) = -iz / e- 2lX( - x - y ) a{y)dy + z 2 / e^^VyM^m^i, A, z)dtdy(3.14) 

J — oo J —oo J — oo 

/■x ry 

m R4 (x,X,z) = 1 + z 2 / / e 2lX{y - t) a(y)a(t)m m (t,X,z)dtdy. (3.15) 



oo *-' — oo 



Iterating the Volterra equations (|3.8j) - (|3.11j) . we prove easily the following lemma 

Lemma 3.1. (i) Setm° L1 (x,X) = 1 and m n L1 {x , A) = ^ x ^ x e 2lX ^~ y ^a{y)a{t)m n L1 1 {t,X)dtdy. Then 
we get by induction 

\m n L1 (x,X)\<— [j a{y)dy) . 

oo 

For x, X € K fixed, the serie m^iix, A, z) = m2i(£, X)z 2n converges normally on each compact subset 

71=0 

of C and satisfies the estimate 



+ oo 



\mLi(x, A, z)| < cosh ( \z\ / a(s)dsV Viet, z e C. 
Moreover, the application z — > m^\{x. A, z) is entire and even. 

(it) Set m° L2 (x, A) = -i e^^^a^dy and m n L2 {x, A) = e ^K^-y) a{y)a{t)m n L2 1 (t, X)dtdy. 

Then we get by induction 

I / f + °° \2n+l 

\mi 2 (x,x)\< (2n+1)i (y • 

oo 

For x, A € K fixed, the serie mL2(x, A, z) = ^""^ ™L2( a; i A)z 2 " +1 converges normally on each compact 

n=Q 

subset of C and satisfies the estimate 

\niL2(%f A, z)| < sinh ( |z| / a(s)ds), Vx e M, z e C. 
Moreover, the application z — > m^iix, X, z) zs entire and odd. 

(tit) Setm° L3 (x,X) =if+°°e 2iX (y-^a(y)dy andml 3 (x,X) = J^°° e- 2 ^ x ~ y )a(y)a(t)m n L - 1 (t,X)dtdy. 
Then we get by induction 

1 / r + °° X 2n+1 

oo 

For x, X € K fixed, the serie mLs{x, A, z) = ^""^ m, L3( x > A)z 2 " +1 converges normally on each compact 
subset of C and satisfies the estimate 

\m L3 (x, A, z)| < sinh ( \z\ f a(s)ds\ \fx 6 R, z e C. 

J X 
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Moreover, the application z — > rriL^x. A, z) is entire and odd. 

(iv) Set m° L4 (x,X) = 1 and m n Li {x,X) = f+°° e- 2l ^ t -y)a{y)a{t)m n L - 1 {t,X)dtdy. Then we get by 



induction 



\mUx,X)\ < 7^_(l +OC a(y)dy) 



For i, Aei fixed, the serie m^ix, A, z) — m^^ix, X)z 2n converges normally on each compact subset 

71=0 

of C and satisfies the estimate 

r+oo 



{itil^x, A, z)\ < cosh ^|z| / a(s)ds S j 1 Vx <E 



z £ 



Moreover, the application z — > m^^x, A, z) is entire and even, 
(v) Note at last the obvious symmetries 



m,Li(x, A, z) — mLi(x, A, z), Vz £ C, 
m L2 (x, A, z) — m L3 (x, A, z), Vz £ C. 

Of course we have similar results for the Faddeev functions m R j(x, A, z), j = 1, ..,4. 
Remark 3.2. Using the usual notations 



(3.16) 
(3.17) 



F L (z,A,z) 



/li(,x,A,z) f L2 {x,X,z) 
f L3 (x,X,z) f L4 (x,X,z) 



, Fr(x, A,z) 



fm(x,X,z) f R2 (x,X,z) 
fm(x,X,z) fm(x, X, z) 



and \3.2\) . we have 



f L j{x,X,z) = e l x m Lj (x,X,z), f Rj (x,X,z) 



i\x 



m Rj (x,X,z), j = l,3, 



f L j{x, X,z) = e 1 x m Lj {x, X,z), f Rj (x, A, z) = e lXx m Rj (x, A, z), j = 2,4 



(3.18) 
(3.19) 



Hence the components fi,j(x,X,z) and f R j(x, A, z) share the same analytic properties as the Faddeev 
functions mi,j(x, A, z) and m R j(x, A, z) respectively. For later use, we mention two additional properties 
of the fLj(x, A, z). First, using (flf. Prop. 2.2) and the analytic continuation, we have 



det{F L (x, A, z)) = 1, Vi€t, z £ C. 



(3.20) 



Second we notice that the fLj(x,X,z) and f R j(x, A, z) satisfy second order differential equations with 
complex potentials. Precisely, the components fLj(x,\,z) and f R j{x, A, z), j — 1,2 satisfy 



d 2 a'(x) d 2 2 a'(x)-\ 
dor a{x) ax a(x) J 

whereas the components fhj{x, A, z) and f R j(x, A, z), j — 3, 4 satisfy 

" T~2 + "TTT" + z a i x ) + lA_ TT f = X f- 

dx^ a(x) ax ayx) i 

The differential equations \3.21\) and i3. 22]) follow directly from the uncoupled integral equations i3.8\ 



(3.21) 



(3.22) 
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Let us now extend the previous result to the matrix Al(\, z). From (|2.21[) and (|3.5|) . we see first that 
the components of Al(\, z) can be expressed by means of the Faddeev functions mLj{x, X, z) as 



ai,i (A, z) = 1 — iz a(x)mL3(x, A, z)dx, (3.23) 
Jr 

a_L2(A,z) = — iz / e~ 2lXx a(x)m,LA(x, A, z)dx, (3.24) 
Jr 

a L3 (\,z) = iz e 2lXx a(x)m L i(x,X 1 z)dx, (3.25) 
Jr 

ol4(A, z) = 1 + iz I a(x)mL2(x, A, z)dx. (3.26) 



Hence we get using Lemma 13.1 1 

Lemma 3.3. (i) For AgI fixed and all z € C, 



x|z 2 "+ 2 , 



an(X,z) = 1 — 1 1 2^ ( / a{x)m r l 3 {x,X)d: 
„=o V - /r 

oo - 

a L 2{\z) = / e- 2iXx a{x)m n L4 (x,X)dx)z 



n=0 



2n+l 



oo „ 

a L3 (X,z) = iV ( / e 2iXx a(x)ml 1 (x,X)dx)z 2n+1 
n=o V « 7 

oo „ 

ciL4(A, z) = 1 + ( / a(x)m r l2(x, X)d. 



x|z 2n+2 . 



(zi^ Se£ A = j a(x)dx. Then 



\a L1 (X,z)\, \a L4 (X,z)\<cosh{A\z\), VzeC, (3.27) 

\a L2 (X,z)\, \a L3 (X,z)\ < wah(A\z\), VzeC. (3.28) 

(Hi) The functions Of,i(A, z) and aLi(X, z) are entire and even in z whereas the functions aL2(X,z) and 
»L3(A, z) are entire and odd in z. Moreover they satisfy the symmetries 



a,Li(X,z) = ciL4,(X,z), VzeC, (3.29) 
a L2 (A,z) = a L3 (X,z), VzeC. (3.30) 

(iv) The following relations hold for all z e C 



a L i(X,z)a L1 (X,z)-a L3 (X,z)a L3 (X,z) = 1, (3.31) 
a Li (X, z)a i4 (A, z) - a L2 {X, z)a L2 (X, z) = 1. (3.32) 

Proof. The properties (i)-(iii) are direct consequences of Lemma [3.11 and of formulae (|3.23[) - (|3.26[) . Let 
us prove (iv). From (|2.23[) . recall that we have for all zel. 

|a L1 (A,z)| 2 -|a L3 (A,z)| 2 -l. 
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This last equality can be rewritten as 



Oli(A, z)a L1 (X, z) - a L3 (X, z)a L i(\, z) = 1, Vz e K. (3.33) 



Since the the functions ai,j(A, z) and a,Lj{\, z) are entire in z by (iii), the equality (|3.33[) extends an- 
alytically to the whole complex plane C. This proves (|3.3ip . Using (|3.29[) . (|3.30l) and (I3.3ip . we get 
(|332l) . " □ 

At this stage, we have proved that the components of the matrix Al(\, z) are entire functions of 
exponential type in the variable z. Precisely, from (|3.27[) and (|3.28l) . we have 

\a Lj (\,z)\<e A W, VzGC, j = l,..,4, (3.34) 

where A = / a(x)dx. We now use the relations (|3.31|) . (|3.32p and the parity properties of the ai,j(A, z) 
to improve this estimate. 

Lemma 3.4. Let A G R be fixed. Then for all z G C 

|azy(A,z)| < e^WI, i = l,..,4. (3.35) 

Proof. From p.3ip . we have 



a L i(A, z)a L i(A, z) - a L3 (A, z)a L3 (A, z) = 1, Vz G 
In particular, we get for purely imaginary z = iy 



oli(A, iy)a L i{\, -iy) - aL3(A, iy)a i3 (A, -iy) = 1, Vy G M. 

But the parity of oli and the imparity of clls yield 

\a L1 {\iy)\ 2 + \aU\ 1 iy)\ 2 = 1, Vy G M. (3.36) 

Similarly, using (I3.32[) instead of (|3.3ip . we get 

\a L2 (\,iy)\ 2 + \a Li (\,iy)\ 2 = 1, Vy G M. (3.37) 

Hence we conclude from (|3.36[) and (13.37)) that 

|0Li(A,iy)|<l, VyeR, j = l,..,4. (3.38) 

Now the estimate (|3.35p is a direct consequence of the Phragmen-Lindclof theorem (see [7], Thm 1.4.3.) 
together with (|3.34[) and p.38p as well as the parity properties of the aLj{\ z). □ 

For later use, we mention that we have the corresponding estimates for the Jost functions fhj{x, A, z) 
and fnj(x, A, z). Precisely 

Lemma 3.5. For all j — 1, .., 4 and for all 

\f Lj (x,X,z)\ < Ce^WIt^)^, (3.39) 
\f Rj (x,X,z)\ < Ce 1 ^ 1 ^^'*. (3.40) 
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Proof. Using the relations (|2?19|) and (|2l20]> . we can prove as in the preceding lemma that the Jost 
functions are bounded on the imaginary axis z G iM, i.e. 

\f Lj (x,X,iy)\, \f Rj (x,X,iy)\ < 1, Vy € R. 

We then conclude by the same argument as above. □ 

We end up this section studying the special case A = in which explicit calculations can be made. 
We have 

Lemma 3.6. Set X = J_ a(s)ds and A = J R a(x)dx. Then 

, . , / cosh(zX) — isinh(zX) \ w _ ,„ 
M R (x,0,z)=\ . . A „( , VzeC, (3.41) 

K ' ' ' \ ismh(zX) cosh(zX) J ' v ' 



/ cosh«A-X)) -isinh(z(A-X)) \ 
M ' U ' J ~ ^ isinh(z(A-X)) cosh(z(A - X)) y ' V 



zeC. (3.42) 



. , , cosh(Az) — isinh(Az) , , , , 



Moreover 

i sinh(Az) cosh(Az) 

Proof. We just prove a few equalities, for instance mL 3 {x, 0,z) = — isinh(z(A — X)) and az,i(0,z) = 
cosh(zA). From Lemmata 13. II and 13.31 recall that 



m L3 (x,0,z) = ^m2 3 (x,0)z 2 " +1 , 

n=0 

a>Li(0,z) = l-i*y]( I a(x)ml 3 (x,0)d. 



n=0 



x)z 2n+ \ 



where 

J m£ 3 (x,0) = ij x °°a{y)dy, 



m n L3 (x,0) = C°° f+°° a(y)a(t)mr 3 \t,0)dtdy 

i / f + °° \2n+l 

ml 3 (x,0)= (2n + 1 J <a)ds) , VneN. 



By induction, we see that 

(2n + l)l\J x 
Therefore, on one hand, we obtain 

oo 

m L3 (x, 0, z) = ^2 m L3 (x 7 0)z 2n+1 = i sinh(z(^l - X)), 

71=0 

and on the other hand, we get 

oo 1 

a L1 (0, *) = ! + £ 7^2Y A2n+2z2n+2 = cosh(Az). 
The other equalities in (|3.41[) . (|3.42[) and (|3.43[) are obtained similarly. □ 
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Remark 3.7. It follows from Lemma \3.6\ that the partial scattering matrices at the energy A = are 
given explicitely by 



cosh(nA) 



and thus depend only on the parameter A. As a consequence, we conclude that the full scattering matrix 
5(0) does not determine uniquely the three parameters of the black hole. 

3.2 Nevanlinna class and uniqueness results 

In this section, we prove that the coefficients aij(A, z) belong to the Nevanlinna class when restricted to 
the half plane II + = {z G C : Re{z) > 0}. As an application, we prove the uniqueness results mentioned 
in the introduction. 

Recall first that the Nevanlinna class iV(II + ) is defined as the set of all analytic functions f(z) on II + 
that satisfy the estimate 

,l_ re i*>. 



sup / kr 

0<r<l , 



1 + re l v 



dip < oo, 



where In (x) — I „ ' . — ' We shall use the following result implicit in l34l . 

[ 0, lux < 0. 1 1 

Lemma 3.8. Let h £ H be an holomorphic function in H + satisfying 

\h{z)\<Ce ARe{z \ Vzen+, (3.44) 
where A and C are two constants. Then h G iV(II + ). 

Proof. Without loss of generality, we can always assume C > 1. Thus using (|3.44j> and ln + (a6) < 
ln + (a) + ln + (6) for a, b > 0, we have 



1^ 



1 - re lv 
1 + re i( f> 



< lnC+ A Re 
lnC*+ A 



1 - re lv 
1 + re^ 
l-r 2 



1 + r 2 + 2r cos <p 

r i - r 2 

Now the well known formula / t. = 27r yields the result. □ 

J_ rr 1 + r z + 2rcosip 

As a direct consequence of Lemmata 13.41 and 13.81 we thus get 

Corollary 3.9. For each X 6 K fixed, the applications z — > a£j(A, z)|n+ belong to N(H + ). 

Let us recall now a usefull uniqueness theorem involving functions in the Nevanlinna class A r (II + ) 
as stated in [34]. We also refer to [38], Thm 15.23, for a conformly equivalent version of this theorem 
involving functions in the Nevanlinna class N(U) where U is the unit disc. 

Theorem 3.10 ([33], Thm 1.3). Let h £ N(U+) satisfying h(n) =0 for all n £ C where C C N* with 

E — = oo. Then h = in II + . 
n 



n 



Hence we deduce from Corollary 13.91 and Theorem 13.101 
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Corollary 3.11. Consider two dS-RN black holes and denote by a^j and a^j the corresponding scattering 
data. Let £cH* satisfying ~S~^ — = oo. Assume that one of the following equality hold 



Then 



a Lj (X,n) = a Lj (X,n), Vn € £, j = l,..,4. 



a Lj (X, z) = a Lj (X, z), Vz G C, j = l,..,4. 



Proof. Apply Theorem 13.101 to the difference functions a,Lj(X,z) — aLj{X,z) which belongs to 7V(n + ) 
thanks to Lemma □ 



In other words, the scattering data a,Lj(X, z) are uniquely determined as functions of z G C from their 
values on the integers n G C. We now extend this uniqueness result assuming that only the reflection 
coefficients L(X,n) or R(X,n) are known as stated in Thm ll.ll Precisely, we prove 

Proposition 3.12. Consider two dS-RN black holes and denote by Z and Z all the corresponding scat- 
tering data. Let £cf satisfying \ — = oo. Assume that there exists a constant c such that one of 

71. 



the following equality hold 



Then 



nec 

L(X,n) 
R(X,n) 



a L1 (X,z) a L2 {X,z) 
a L3 (A, z) a i4 (A, z) 



e- 2lXc L(X,n), 
e 2lXc R(X,n), 



a L1 (X,z) 
^ Xc ~a L3 {X,z) 



Vn G C, 
Vn G C, 

e 2lXc a L2 (X,z) 

&L4(A, z) 



Vz G C. 



Proof. Assume (|3.45|) . Then 

a L3 (X,n)a L1 (X,n) 



-2i\c 



a L1 (X,n)a L z{X,n), Vn G £. 



(3.45) 
(3.46) 

(3.47) 
(3.48) 



By Lemmata 13.41 and 13.81 the product functions a,£ 3 (A, z)a£i(A, z) and a^i(A, z)a£ 3 (A, z) belong to the 
Nevanlinna class N, the equality p.48[) extends analytically to the whole complex plane C. Hence 



Oi3(A, z)a L i(X, z) = e 2lXc a L i(X, z)a i3 (A, z), Vz G 
But recall from Lemma [3~3l (iv), that 



(3.49) 



a/,1 (A, z)a L i(X, z) = 1 - a L3 (X, z)a i3 (A, z), Vz G C, 

from which we deduce that 0£,i(A, z) and ci£ 3 (A, z) have no common zeros. Hence we infer from (|3.49p 
that the zeros of a^j(X, z) and a™ (A, z) for j — 1,3 coincide with the same multiplicity. 

Now recall that the function ai,\[X,z) is even. Thus we can write ati(A, z) = g(z 2 ) where g is an 
entire function. Since cili(X,z) is of order 1 (i.e. \oli(X, z) < e A ' 2 '), we deduce that g is of order |. 
Hence the Hadamard's factorization theorem, (see [7], Th 2.7.1), yields 



ff(o=<?n(i-£), 



n=l 



21 



where the £„ 7^ are the zeros of g counted according to multiplicity, G — g(0) — flLi(A, 0) = 1 by 
Lemma T3.31 (i). But note that ( n = z 2 where the z n are the zeros of Oli(A, z) by definition of g. Hence 
we obtain 

00 2 
a L1 {X,z) = Yl (l- - 
n=l 



Similarly, we have 



1 (A^)=n( i -?2 



where the z n are the zeros of a 1,1 (A, z). Since z n — z n by the previous discussion, we conclude that 

ati (A, z) = ctLi (A, 2), VzeC. (3.50) 

From (|3~50)) . (|3~49)) . (|3~29)) and (fOP]) . we thus deduce (pT47j) . The proof starting from (f3T46|) is analogous 
and so we omit it. □ 

In the same way, if we assume that only the transmission coefficient T(A, n) are known, we have the 
following result : 

Proposition 3.13. Consider two dS-RN black holes and denote by Z and Z all the corresponding scat- 
tering data. Let C C N* satisfying y — = 00. Assume that T(X,n) = T(X,n) for all n G L. Then, 

nec 

there exists a constant c € R such that 



a L x(X,z) a L2 (X,z) 
a L3 (A, z) a i4 (A, z) 



a L1 {X,z) e 2lXc h L2 {X,z) 
a L3 {X,z) a L4 {X,z) 



-2i\c~ 



Vz e C. 



Proof. By Lemma T3.31 (iii) and Corollary 13. 11[ we have 

a Lj (X,z) = a Lj (X,z), Vz £ C j = 1,4. 



(3.51) 



(3.52) 



Now, we use the same strategy as in Prop. I3.12l and we set f(z) = L3 ^ - — -. Using that aL3(X, 0) = 0, 

we see that f(z) is an even entire function of order 1. Thus, we can write f(z) = g(z 2 ) where g is an 
entire function of order \. Using the Hadamard's factorization theorem, we obtain as previously 



00 2 



where 2m is the multiplicity of 0, G is a constant and the z n are the zeros of / counted according to 
multiplicity. From (|3.52p and Lemma 13.31 (iv), we have 



f(z)f(z) = f(z)f(z), 



where f(z) = L3 ^~ - — -, or equivalently 



ri—1 11 ' n— 1 IL 
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It follows that 

| G |=| G |, TO = TO, z n = z„ 

Recalling that ol?,{X,z) = zf(z) and a~L 3 {X,z) = zf(z), we see that there exists a constant c g R such 
that 

a L3 {X,z) = e~ 2lXc a L3 (\,z), z G C. 

Finally, using Lemma l3.3[ (iii) we obtain the result. □ 

We finish this section with a first application of the previous uniqueness results to the study of an 
inverse scattering problem in which the scattering matrix S(X, n) is supposed to be known on an interval 
of energy only (and not simply at a fixed energy A). We shall obtain in an elementary way one of the 
results in [14J . Precisely we prove 

Corollary 3.14. Assume that one of the scattering coefficients T(A,n), L{X,n) or R(X,n) be known (in 

the sense of Theorem for all n G C with } — — oo and on a (possibly small) interval of energy 

neC 

X € I. Then the parameters M,Q 2 ,A of the black hole are uniquely determined. 

Proof. Consider two dS-RN black holes with parameters M, Q, A and M, Q, A respectively such that for 
all n G C and for all X £ I, one of the following equality holds 

L(X, n) = e- 2lXc L(X,n), 
R(X,n) = e 2lXc R(X,n), 
T(X,n) = T(A,n). 

Then, applying Propositions 13.121 and 13.131 we see in particular that 

a L2 (A, z) = e 2tXc a L2 (X, z), Vz e C, VA e I. 

From the first term of the series defining ollz(X, z) and 5^2 (A, z) (see Lemma l3~3l (ii)), we thus obtain 

a(2A) = e 2lXc a{2X), VA G I, (3.53) 

where a and a denote the Fourier transforms of the potentials a and a. Since these potentials are 
exponentially decreasing at both horizons, their Fourier transforms a and a are analytic in a small strip 
around the real axis, i.e. on K = {X G C, |7ro(A)| < e} for e small enough. Thus the equality (|3.53l) 
extends analytically to the whole strip K . In particular, we have 

a(2A) = e 2iAc a(2A), VA G K, 

and therefore 

a(x) = a(x - c), Vz G R. (3.54) 

According to the proof following (I5.18[) in section 5, the equality (|3.54[) is the essential ingredient to 
conclude to the uniqueness of the black holes' parameters. □ 
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4 Large z asymptotics of the scattering data 



In this section, we obtain the asymptotic expansion of the scattering data when the coupling constant 
z — > +00, z real. We emphasize that this is the only place in this paper where we use the exponential decay 
of the potential a(x) at both horizons and thus the asymptotically hyperbolic nature of the geometry. 
The main tool to obtain these asymptotics is a simple change of variable X = g(x), called the Liouville 
transformation. As we shall see in the next section, this Liouville transformation as well as the asymptotics 
of the scattering data will be useful to study the inverse scattering problem. 

As a by-product, we shall obtain a very simple reconstruction formula for the surface gravities n± from 
the scattering reflexion coefficients L(X, n) and R(X, n), n £ N. Recall that the surface gravities have an 
important physical meaning through the celebrated Hawking effect which asserts that static observers 
located far from both horizons measure at late times an isotropic background of thermal radiations 
emitted from the horizons. The rate of radiation is given by the "temperatures" of the cosmological and 
event horizons which are shown to be proportional to the surface gravities n±. 



4.1 The Liouville transformation. 

We follow the strategy adopted by K. Chadan, R. Kobayashi and M. Musette (see [10] . QT|). Considering 
the differential equations fl3.21[) and (|3.22p satisfied by the Jost functions fLj(x, A, z) and fnj(x, A, z), we 
use a Liouville transformation, i.e. a simple change of variable X = g{x), that transforms the equations 
()3.2ip and (|3.22|) into singular Sturm-Liouville differential equations in which the coupling constant z 
becomes the spectral parameter (see Lemma 14.11 below) . 

Let us define precisely this Liouville transformation. We denote 

X = g{x) = [ a{t) dt. (4.1) 
J — 00 

Clearly, g : K — >]0, A[ is a C 1 -diffeomorphism where 

/+00 
a(t) dt. (4.2) 
-00 

For the sake of simplicity, we denote h = g~ x the inverse diffeomorphism of g and we use the notation 
df 

f(X) = -7^{X). We also define for j = 1,...,4, and for X e]0,A[, 

f j (X,X,z) = f Lj (h(X),X,z). (4.3) 

gj (X,X,z) = f Rj (h(X),X,z). (4.4) 

We begin with an elementary lemma which states that, in the variable X, the Jost solutions fj(X, A, z) 
and gj(X, A, z) satisfy Sturm-Liouville equations with potentials having quadratic singularities at the 
boundaries. 

Lemma 4.1. 

1. For j = 1, 2, fj{X, A, z) and gj(X, A, z) satisfy on ]0, A[ the Sturm-Liouville equation 

y" + q(X)y = z 2 y. (4.5) 
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2. For j = 3, A, fj{X, A, z) and gj(X, A, z) satisfy on ]0, A[ the Sturm- Liouville equation 

y" + qJX)y = z 2 Vl (4.6) 
A 2 a' (a;) 

where the potential q(X) = X 2 h'(X) 2 — iXh" (X) — + iX has the asymptotics 

a (x) a a [x) 

q(X) = —r+q-iX), with w_ = ^ + i— , and q_(X) = 0(1), 0. (4.7) 



y(.V] ( A U * X )2 + q+i - X) ' With W+ = ^ +i ZT i, " <1 <U ' X,: X ""' A ' ' 4 ' 



Proof. Using A3. 21ft . (13.221) and = a(x), a straightforward calculation gives (14.5ft and (14.61) . Let us 

ax 

establish (14.71) . By Lemma T2. 11 recall that the potential a has the asymptotics when x — > — oo 

o(x) = a_ e*"* + 0(e 3K - x ), 
a'(x) = K_a_ e K -" + 0(e 3K - x ). 

When i ^ -co, or equivalently when X — »• 0, it thus follows that 

/X -t 
a (t) dt= — e K ~ x + 0{e 3K - x ) = — a(x) + 0(e 3n - x ), 
-oo K - K - 

from which we also have 

e K - x = 0(X), X^O. 

Hence we find 

a(x) = k^X + 0(X 3 ), X^O. (4.9) 

Simarly we get 

a'(x) = k 2 _X + 0(X 3 ). X -> 0. (4.10) 
Now, (l4~9l and (l4T0|) yield immediately (j4~7) . The proof of (|4~8)l is similar. □ 

In the next lemma, we calculate the wronskians of some pairs of Jost functions in the variable X. We 
recall that the wronskian of two functions /, g is given by W(f, g) = fg' — f'g. These wronskians will be 
useful to obtain the asymptotics of the scattering data and to solve the inverse problem. 

Lemma 4.2. For z 6 C, we have : 

W(hJ 2 ) = W( gi ,g 2 ) = W(.hJ 4 ) = W(g 3 ,g 4 )=tz. 

Proof. For example, let us calculate W{f\,f2). Using —j— = a(x) again, it is clear that 

W(hJ 2 ) = -}-W(f L1 (x,X,z), f L 2(x,X,z)). (4.11) 
a{x) 

Using (|3.ip and (|3.20[) . we obtain easily : 

W{h,h) = iz(f L ifLi - /w/w) = iz det F L = iz. 

□ 
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4.2 Asymptotics of the Jost functions. 

Singular Sturm Liouville equations like (|4.5[) - (|4.8|) have been yet studied by Freiling and Yurko in 
[18) . In particular, fundamental system of solutions with precise large z asymptotics have been obtained 
there. In theory, we could thus establish a dictionary between our Jost functions and these fundamental 
solutions to infer the asymptotics of fj(X) and gj(X). However, the asymptotics of the fundamental 
solutions of |18j are given up to multiplicative constants that, as regards the inverse problem we have in 
mind, require to be determined precisely. Since this is not a straightforward task, we prefer to follow a 
self-contained and elementary approach that is similar to the usual perturbative argument for the Sturm 
Liouville equation (|4.5lFI . but only uses the series expansion for the Jost functions obtained in Section [3l 
Note that this is not at all compulsory but allows us to keep track of the different multiplicative constants 
appearing in the perturbative method from the very begining. 

Let us study first the Jost function fzi(x, A, z). From Lemma l37[l and (|3.18p . we recall that for z e C, 

+00 

f L1 (x,X,z) = e iXx Y, m li( x > X ) z2n ( 4 - 13 ) 

n=0 

r+00 r+00 

m2i(a;,A) = / / e 2lX{t ' v) a{y) a{t) m^t.A) dt dy , for n > 1. 

iateh 

+00 



where 
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So, using the Liouville transformation, we immediately obtain 
Lemma 4.3. 

fi(X, X,z) = Y a n(X, A) z 2n , (4.14) 



n=0 

where 

ao(X,X) = e^W, 

a n (X,X) = e iXh W f£fA e -2i\h(Y) e iXh(T) an _^T,X) dT dY , for n > 1. 1 ' 

We shall consider fi(X,X,z) as a perturbation of a function fx(X, X, z) where fi(X,X,z) is given 
by the same serie (|4.14j) - (|4.15j) as fi(X, A, z) with h(X) = g^ 1 (X) replaced by another diffeomorphism 
denoted h + (X) = g^_ 1 (X). The diffeomorphism g+ in turn is defined in the same manner as g but we 

replace a(x) by its equivalent at +00. More precisely, if we write A — g(x) = J + °° a(t) dt, it is natural to 
set for X e]0,A[ 

r+00 „ 

A - g+(x) = / a+ e K+t dt = ±e K + x . 

Jx « + 

So, we define 

h+(X) = g+\X) = — ]og(A -X)+ C+, (4.16) 

K+ 



3 A fundamental system of solutions (FSS) of | |4. 511 - 114. 8H could be constructed by perturbation of solutions for the free 
equation 

y" + ^y = z2 v- ( 4 - 12 ) 

Noticing that the modified Bessel functions form a FSS for 1 14. 121 . we could thus construct a FSS for the full equation using 
good estimates on the associated Green kernel (see |18i 1251 137] where a similar procedure can be found). In our approach, 
we shall retrieve directly that the Jost functions are perturbations of the modified Bessel functions (see Proposition 14. 12l 
and thus deduce their asymptotics from the well known asymptotics of the lattcrs. 



2G 



with 

C+ = — log(-^). (4.17) 

f+ 



The function (X, A, z) is thus given by 



+00 



f+(X,\,z) = Y,a+(X,X) z 2n , (4.18) 



n=0 

where 



a+(X,A) = e^ + W 

a+(X,A) = e 4A,1+ W / A ( A e -^x h +(Y) e iXh+(T) a +_ i(T , A) dT dr , for n > 1. 
Jx Jy 

Thanks to our choice of diffeomorphism h + , the coefficients of the serie f|4. 18[) can be explicitely 
calculated. Precisely, denoting by T the well-known Gamma function, we have 

Lemma 4.4. For X e]0, A[, z e C and for all n > 

a+(*,A) = (- — )^ r(l-^ + ) , 2 „ r . 1 r— (A-Xf n+ ^ 

a + 2 zn 1 (n + 1 — 1/+) n! 

^ + = l-zA (4.19) 

Proof. We prove the formula by induction. For n = 0, the result is clear by (|4.16l) . f|4. 1T|) . For n > 1, an 
elementary calculation gives 

a+(X, A) = (-— )^ =~ =~ (A - X) 2n+i T (4.20) 

nV ' ' a+> (1 + 2»A)...(2n-l + 2sA) 2---(2n) V ; V ; 



Using the functional equality T(z + 1) = zT(z), Lemma \A. 41 is proved. □ 

Now, it turns out that the expressions for the coefficients a+(X, A) can be written in terms of the 
modified Bessel function I_ u (x). Let us recall its definition ([H]) Eq. (5.7.1), p. 108), 

7 -^) = E vl — 1 (?) ' ^ € C, | Ar 9 x |< tt. (4.21) 

' l(n— v + lln V2/ 

n=0 V ; 



We deduce 
Corollary 4.5. 

L For X 6]0,A[ and z g 



f+(x,x,z) = (-^±)^ r(i - Va^T (£)-+ j_„ + (z(^-a:)). (4.22) 
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2. Let X\ e]0, A[ fixed. Then, for k = 0,1 and for all X <=.]§, X\[, the following asymptotics hold when 
z — > +00, z real: 

f^ k \x,\,z) = (-!)* BT1 (_*t)£ r(l-* + ) (l + 0(-)). (4.23) 

VZ7T O-l- v 2 / 



X 

Proof: The first assertion comes from Lemma (|4.4j) and (I4.21j) . observing that (— )" I-v{ x ) is holo- 

morphic on C. For the second one, let us recall the well-known asymptotics for the modified Bcsscl 
function I- V (x), v G C, k = 0, 1, when x — > +00 : 

/^) = i (1 + 0(|*|- 1 )) , (4.24) 
V 27ra; 

For the case A; = 0, we refer to ([25J, Eq. (5.11.8), p. 123). The case k = 1 follows from the previous case 
together with the identity (see [25], Eq. (5.7.9), p. 110) : 

Then, the asymptotics (|4.23l) are a simple consequence of (I4.24p . ^> 

Remark 4.6. Let us study the special case A = 0. From the definitions, it is immediate to see that 
fi(X,0,z) — fi(X,0,z) and v+ = 5. So, Corollary \4-5\ entails that 



h{X,\z) = T{\) \J Z[A 2 X) I_ h {z{A-X)). (4.25) 



Using that r( — ) = \fn and the equality (see f29f , Eq. (5.8.5), p. 112), 



I-i(x) = \ — coshx , (4.26) 

2 \ TTX 



we get fi(X, 0, z) = cosh(z(A — X)). Hence we rediscover the result obtained in Lemma \3.6l 

In order to estimate fi(X, A, z) — f^(X, A, z), we have to control h(X) — h + (X). Since in the construc- 
tion of h+(X), we have replaced a(x) by its asymptotic at +00, it is hopeless to get globally estimates 
on ]0, A[. This is why we shall work in ]Xq, A[ where Xq g]0, A[ is fixed. We have the following result : 

Lemma 4.7. Let Xq e]0, A[ fixed. Then, there exists C > such that for k = 0, 1, 2 

I ftW(I)-/)<"(I) I < C{A-Xf- k , Vl€]I ,4 (4-27) 

Proof. We only give the proof for k = since the other cases are similar. For X G]X , A[, by Lemma |2~T1 
we have 

r+00 p+00 

A-X= a(t)dt= (a+e K+ * + 0(e 3K +*)) dt = ~-!-e K + x + 0(e 3K+x ). (4.28) 

Jx Jx K + 



Then, 



log(-^±) + log(A - X) = K+ x + 0(e 2K + x ). 
a + 



By (|L28|) . we see that 0(e 2K+x ) = 0((A - X) 2 ). So, using (14J7) . we obtain : 

x = h{X) = — \og(A -X) + C+ + 0((A - X) 2 ) = h+(X) + 0((A - X) 2 ). 



□ 
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We can deduce from Lemma T4. 71 some useful properties for a n (X, A) and a„(X, A). 
Lemma 4.8. 

1. For n > and X g]0, A[, we have 



(A~X) 2n 
~W) 



an(X,X) |< V ( , (4.29) 



& Forn>0 andIe]Io,i[, 

<(jr,A) = 0((A-X) 2 "- 1 ). (4.30) 

TTie same estimates hold for a+(X, A). 

Proof. The first point is clear by a simple induction. Let us prove the second one. For n = 0, we observe 

that Oopf, A) = ih'(X) ao(X). But, since /i+(X) = 1—. r, Lemma [4771 yields the estimate 

k+ [A — X J 



(4.31) 



which proves the result. For n > 1, we have 



a n (X, A) = iX h'{X) a n (X) - e -' Ah W / e -Wn ajj _ 1 ( T) A ) dT . 



v 



Hence the result follows from (14.29[) and (|4.31[) by induction. For a+(X, A), the proof is identical. □ 
Now, we want to control the difference fi(X, A, z) — ff{X, A, z). To do this, we set 

e n (X,X) = a n (X,X)-a+(X,X), (4.32) 

and thus, we have 

h(X, A, z) - f+(X, \z) = J2 A) z 2n . (4.33) 



In the next lemma, we show that a n (X, A), a+(X, A) and e n (X, A) satisfy second order differential equa- 
tions. 

Lemma 4.9. 

1. For n > 1, a n {X, A) and a£(X, A) satisfy on ]0, A[, 

a';{X,X)+q(X) a n (X,X) = o„_i (4.34) 

a+"(X,A)+ g+ (X) a+(X,A) - ^,(1,)), (4.35) 

w/iere 

= X 2 h + '(Xf iXh + "(X) = + -£) {A ^ X)2 - (4.36) 
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2. For n > 1, e n (X, A) satisfies on ]0, A[, 

e n "{X,X)+q+{X) e n (X,X) = e n . 1 (X,X)-q (X) a n (X,X) (4.37) 

where 

q (X) = q(X) - q + (X). (4.38) 

+00 

Proof. Since h(X,\z) = ^ a n {x,X) z 2n , (|434l) follows directly from g3]). The proof of (|4351) is 

n=0 

identical. At last, (|4.3T[) is a direct consequence of (14.34[) and (|4.35j) . □ 

Now we show that the equation (I4.37[) can be rewritten as an integral equation which will be useful 
to estimate the error term e n . 

Lemma 4.10. 

1. qo(X)eL°°(X Q ,A). 

2. For n > 1 and for all X £]Xo, A[, e n {X, X) satisfies the integral equation 

e n (X, A) = e iXh +W [ A [ A e ~ 2iXh +^ e lXh +^ [e»_i(T, A) - q (T) a n (T, A)] dT dY. (4.39) 
Jx Jy 



Proof. The first point is clear by Lemma 14.71 after having noticed that 

q (X) = X 2 (h'(X) - h + '(X)) (h'(X) + h + '{X)) - iX(h"(X) - h+"(X)) (4.40) 
To prove the second point, we denote by f n (X, A) the (R.H.S) of (|4.39|) . Using (|4.37[) . we have : 

f n (X, X) = e iXh +W f f e -2iXh+(Y) e iXh+(T) gn // (T) A) dT dy 
Jx Jy 

+ e iXh +W [ ( A e -2r\h + (Y) e iXh + (T) ?+(T) 6n(TjA) dT dy _ 

Jx Jy 

Integrating by part twice the first integral and using Lemma T4.8I yield (|4.39|) . □ 

In the next lemma, we estimate e n {X, X) and its derivative. 
Lemma 4.11. There exists a constant C > such that for all n > and for all X £]Xq, A[, 

(A - Y) 2n + 2 

I e n (X, A) | < C(n+ 1) 1 (4.41) 

(A — x) 2n+1 

I e n '(X, A) | < C (n + 1) 1 (4.42) 
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Proof. We prove the lemma by induction. For n — 0, we have 

e (x, A) = e lXh{x) - e iXh+{x) = e lXh+{x) (^KHx)~h + (x)) _ ^ _ 
So, by Lemma I4T7I we obtain for a suitable constant C > 0, 

\e (X,X) |< | (A-X) 2 , VXe(X ,A). 
Assuming that the property is true for n — 1 , we have by (|4.29j) and (|4.39|) , 

e n (X,X)\< [ [ (Cn (A 2 " + c (A ^„T, )2n > ) dT <^ VIe(I 0l A), 



(2n)! (2n) 



where we have supposed that C >\\qo ||i=o(x ,A)- So, 

\e n (X,\)\<C(n + l) { ; , VIe(I ,A). 

We prove (|4.42l) similarly. □ 

Now, we can establish the main result of this section : 
Proposition 4.12. 

1. There exists C > such that for k = 0, 1, for all X &\Xq, A[ and all z > 0, 

| fi k \x,X,z)-f+ ik \x,X,z) | < C {A-X) z k - x e z ^ x \ (4.43) 

2. For fixed Xq < X\ with Xj e]0,A[, Vfc = 0,1, \/X g}Xq,Xi[, we have the following asymptotics, 
when z — > +oo, 

f[ k \x,\,z) = (-l) k 2 ^L (-^±)# r(l-i/+) (l + 0(i)). (4.44) 

V27T a + V z / 

Proof. We only prove (|4.43l) in the case k = since the case A; = 1 is similar. By Lemma 14.111 and for 
z > 0, we have : 

\h(X,\,z)-f+(X,\,z)\ < I e n(x,X) | z 2n 

< g (a x) y (" + 1 ) (^-^)) 2 " +1 

" z ( ' ^ 2n + 2 (2n + l)! 

n— 

< — (A — X) sinh(z(A - X)) 
2z 

< g (A-X)e< A ~ x \ 
z 

Now, since z is real, (I4.44j) follows from Corollary 14.51 and (|4.43l) . □ 
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In order to calculate the asymptotics of f 2 (X, A, z), we follow the same method as for fi(X, A, z). We 
thus only give the final results omitting the details. First, we construct f 2 (X, A, z ) which approximates 
f 2 (X, A, z) as in We obtain : 



f+(X, A, z) = -i (-^) »+ T(l - M+ ) VA^T / x ( Z (A - X)), (4.45) 

z 

where 

M+ = I+jA. (4.46) 

2 K + 

Then, using the well-known asymptotics for the modified Bessel functions, we deduce : 

Proposition 4.13. For fixed Xq < X\ with Xj £]0,A[, Vfc = 0, 1, VX g]Xq,Xi[, we have the following 
asymptotics, when z — > +oo, 



/f(X ) A,.) = (-lf +1 i 2 -Z r(l- M+ ) z k+ ^ e«) (l + 0(-)). (4.47) 

Remark 4.14. As previously, let us study the special case A = 0. We have f 2 (X, 0, z) = f%{X, 0, z) and 
fi+ = \. Hence, using the equality (see [29], Eq. (5.8.5), p. 112) 

Jx(x) = J^- sinha;, (4.48) 



we find that f 2 (X, 0, z) = —i sinh(z(A — X)) as expected according to Lemma \3.b\ 

Remark 4.15. Note that the asymptotics given in \4-44ty an d \4-4ty only depend on n + and a + , i.e. on 
some parameters of the black hole at the cosmological horizon X — A (see \l-4ty and \2. 7| )J. This is natural 
since the Jost functions fi(X, A, z) and f 2 (X,X,z) are solutions of |^.5[ ) with boundaries conditions at 
X = A and since we work in the range X > Xq > 0, i.e. far from the event horizon. We emphasize 
then that the geometry of the black hole at the event horizon X = does not appear in these asymptotics. 
Of course it won't be the case for the scattering data. We also stress the fact that we have only assumed 
that X < X\ < A by convenience: the asymptotics of the derivative of the Jost functions fj(X, A, z) are 
simpler under this condition. 

In order to obtain the asymptotics of the scattering data, we need to calculate the asymptotics of 
the Jost functions gj(X, A, z). Since the procedure is the same as the one for the fj(X,X,z), we give 
without proof the main steps to obtain the asymptotics of gj(X, A, z), j = 1,2, when z — > +oo. Since 
gj(X, A, z) satisfies (|4.5|) with a boundary condition at X = 0, we work with an other diffcomorphism, 
denoted by h-(X), in order to construct the functions g~(X,X,z) that approximate gj(X, A, z). This 
new diffcomorphism is defined as follows : 

h- (X) = — log X + CL , (4.49) 

K- 



C- = — log ( — ) . (4.50) 



where 

C- = 

K~ \ a- 

As previously, we can calculate explicitely gj (X, A, z) and we easily obtain the following equalities 

g -(X,X,z) = (— )- VXT(1-^) C-y- I. v _{zX) , 
a_ 2 
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where 

v- = \-i±, (4.51) 

and 

to(x,\,z) = i v^r(i- M _) { z -y- h-^(zx) , 



where 



fi- = -+i— ■ (4.52) 

Z K — 



Remark 4.16. In the special case X — 0, we have gj(X, 0, z) = (X, 0, z) and \i- = i/_ = |. Using the 
formulae {4-26\j and \4-4<ty f or ^ e modified Bessel functions of half-integral order, we find g^ (X, 0, z) = 
cosh(zX) and g^{X, 0, z) — i sinh(zX) (see Remark\3.~ 



The g± (X, A, z) are perturbations of the gj(X, A, z). Precisely, we have 
Lemma 4.17. For X x e]Q, A[ /tared, t/iere exisis C > such that Vfc = 0, 1, VAT G]0, Xi [, Vz > 0, 

| <?f (X,A,«)- fl 7 (k) (A:,A,z) | < C X z*- 1 e zX (4.53) 
Then, using the asymptotics of the modified Bessel functions (|4.24|l . we obtain : 

Proposition 4.18. For fixed Xq < X\ with Xj G]0,A[, Vfc = 0, 1, VX £]Xq,Xx[, we have the following 
asymptotics, when z — > +oo, 

g [ k \x,\,z) = ^= (— )- r(l-i/_) z*"- e** (l + 0(i)), (4.54) 
V27T a- v z / 

A, z) = i HZl (^)- - r(l - e zX (l + O(-)) , (4.55) 

4.3 Asymptotics of the scattering data. 

In this section, we put together all the previous results and calculate the asymptotics of a^j(X, z), j 

1 . 4 when z — > +oo. First, we recall that for all id, 

F L (x,\,z)=F R (x,\,z) A l (X,z). (4.56) 

Calculating (|4.56|) components by components, it follows that (in the variable X) 

fi(X, A, z) = a L1 (X, z) gi(X, A, z) + a L3 (X, z) g 2 (X, A, z) 

fi{X, A, z) = a L2 (X, z) gi(X, A, z) + a L4 (A, z) g 2 (^, A, z) 
So, by Lemma l4~2l we obtain for z 7^ : 

«li(A,z) = — W(f 1 ,g 2 ) , a L2 (A,z) = — W{f 2 ,g 2 ), 

iz iz 

a L3 (X,z) = -—W(fi,gi), a L4 ,(X,z) = -— W(f 2 ,gi). 

iz iz 

The following theorem is an easy consequence of Propositions 14.121 14.131 and 14.181 : 
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Theorem 4.19. When z — >• +00, we have : 



1 




27 


V a + 


— i 


/ '_K+ 


27 


\ a + 


i 


f_K+ 




V «+ 


1 




27 


V a+ 



K- 



1 «A\„/i u\ rzvHit-it) 



MM- - I--) ■(^)--r(l + ^)r(l-^)(£)«<*-*V(«o, 



We deduce from Theorem 14. 191 the asymptotics of the transmission and reflexion coefficients T(A, 2), 
L(A, 0) and R(X, z). From the definitions (|2.25|) of these coefficients, we get 

Theorem 4.20. When z — » +oo ; we have : 



UJ UJ r (±-£)r(i + £) ^ 
(* + £) 

_ iA \ 

(4 -ft) 



UJ r (|-£) V2 



— r ( i - ^ ) 2iA 



Remark 4.21. ^4s expected, the asymptotic of the transmission coefficient T{X,z) depends on the pa- 
rameters k± and a±, i.e. on the geometries of both event and cosmological horizons. On the other hand, 
the asymptotic of the reflection coefficient L(X,z) depends only on k_ and a_ - the geometry of the event 
horizon - whereas the asymptotic of R(X, z) depends on k + and a + - the geometry of the cosmological 
horizon. 



4.4 Reconstruction formulae for k± 

As a by-product of the asymptotics obtained in Theorem l4.20[ we find simple reconstruction formulae for 
the surface gravities k± from the scattering reflexion coefficients L(X,n) and R(X,n), neN. As already 
mentioned at the beginning of this section, these quantities are meaningful in the Hawking effect. 

Theorem 4.22. For all p E N, we have : 



L(X,pn) 
hm — — — - 

L(X,n) 

hm 

n^+oc R(X, n) 
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We can now determine easily the surface gravities for nonzero energies A ^ 0. For example, if we set : 

-MA logp 

Up = e K - , (4.ol) 

which is known, we obtain when p — > +00 : 

^±i = 1-^ + 0(1). (4.62) 

U p K-P p z 

This permits to calculate k_ if A 7^ 0. 

5 The inverse scattering problem. 

In this section, we prove our main result Thm fTTTl that is we prove the uniqueness of the potential a(x) 
up to translations as well as the uniqueness of the parameters (M, Q 2 , A) of a dS-RN black hole from the 
knowledge of either the transmission coefficient T(X,n), or the reflection coefficients L(X,n) or R(X,n) 
at a fixed energy X ^ and for all n £ C C N* satisfying the Miintz condition J2 n ec h ~ 00 ■ 

Consider thus two dS-RN black holes with parameters (M, Q, A) and (M, Q, A) respectively. We shall 
denote by a(x) and a(x) the corresponding potentials appearing in the Dirac equation and satisfying the 
hypotheses of section 2. We shall also use the notation Z for all the scattering data associated with the 
potential a. As explained in the introduction, we assume that there exists a constant c such that one of 
the following equalities hold for all n g /j3 

T(X,n) = f{X,n), 

L(X,n) = e" 2iAc L(X,n), (5.1) 
R(X,n) = e 2lXc R(X,n). 

By Propositions 13.121 and 13.131 we deduce from (|5.ip that 

oli(A, z) = a L1 (X, z) , a L2 (A, z) = e 2lAc 5 L2 (A, z), (5.2) 
a L3 (X,z) = e~ 2lXc a L3 (A, z) , a L4 (A, z) = a L4 (A, z). (5.3) 

Thus, it follows from the asymptotics of Theorem 14. 191 that : 

/+00 r + co 

a{x) dx = a(x) dx = A. (5.4) 

-00 J — 00 

Hence, we can define the diffeomorphisms h, h : ]0, A[— > R as the inverses of the Liouville transforms g 
and g given by (|4.ip in which we use the potentials a(x) and a(a;) respectively. 

Now, following a strategy relatively close to jTB], we introduce, for X e]0, ^4[, the matrix 



P(X,X,z) = 

defined by 



Pi{X,X,z) P 2 (X,X,z) 
P 3 {X,X,z) Pi{X,X,z) 



P(X, A, z) F R (h(X), A, z) = F R (h(X), A, z) e lXcTl , (5.5) 



4 Recall that we add a constant c in 115.11 1 to include the possiblity of describing the same dS-RN black hole by two 
different RW variables and make our result coordinates independent. 
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where Fr — and Fr — are the Jost solutions from the right associated with a(x) and a(x). 

To simplify the notation, for k = 1, ...,4, we set as in Section 4: 

f k (X, A, z) = f Lk (h(X), A, z), / fe (X, A, z) = f L k(h(X), A, z), 
fffc(X, A, z) = f Rk (h(X), A, z), fa(X, A, z) = f Rk (h(X), A, z). 

Using that det Fr = 1 and det Fr = 1, we obtain the following equalities : 

' Pi(X,X,z) = e lAc g x g 4 ~ e- lXc g 2 g 3 , 
P 2 (X,X,z) = -e lXc 9l ~g 2 -e- lXc g 2 gi- [ 1 

It follows from (|5.6[) and the analytical properties of the Jost functions that, for j = 1,2, the applica- 
tions z —¥ Pj(X, A, z) are holomorphic on C and of exponential type. Moreover, by Lemma 13.51 these 
applications are bounded on the imaginary axis iWL. 

We shall now prove that the applications z — > Pj(X, A, z) are also bounded on the real axis. To do 
this, we first make some elementary algebraic transformations on Pj(X, A, z). We write : 

P X (X, A, z) = e* Ac 91 g 4 + e* Ac 9l (g 4 - g 4 ) - e~ iXc g 2 g 3 - e' lAc g 2 (g 3 - g 3 ), 
= e jAc ( 91 g 4 - g 2 g 3 ) - e~ lXc g 2 (g 3 - e 2lAc g 3 ) 

+ e lXc 9l (34 - .94) - e- tXc 52 (33 - 9a), 
= e lXc + e iAc 32 33 - e^ Ac 32 33 + e lAc g x (34 - 34), 

where we have used that det Fr = 31 34 — 32 33 = 1. Since Fl{x, A, z) = Fr(x, A, z) Aj,(A, z), we get 
using (|5.3j) : 

.94 = (/* - «L2 33), 34 = (/4 - e~ 2lXc a L2 33). 

So, we obtain immediately : 

P X (X, A, z) - e 4Ac + e lAc (/ 4 - / 4 ) + (32 + a L2 ) (e lAc 33 - e" lAc §3). 

«L4 CIL4 

Using again Fl(x, A, z) = Fr(x, A, z) Al(A, z), we see that f 2 = a^ 2 31 + g 2 . Thus, we get 

P X (X, A, z) = e lAc + e* Ac (/ 4 - / 4 ) + — (e lAc 33 - e-* Ac .93). (5.7) 

»L4 CIL4 

Similarly, /^(AT, A, z) can be expressed as : 

P 2 (X, A, z) = — (e- jAc / 2 .91 - e lAc / 2 31) • (5.8) 

OL4 V / 

We shall now use some estimates obtained in the previous sections. First, it follows from Lemma l3~5l 
that for z > and for all j = 1, .. , 4 : 

\fj(X,X,z)\, \MX,X,z)\<e^ A ^ , Iff^X.A^)!, | A, z) |< e zX . (5.9) 

Second, using Theorem 14.191 it is easy to see that for z real and large enough 

\a L i(X,z)\>Ce Az , z»l. (5.10) 
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Hence, using (|5.7|) . ()5.8j) . (|5.9j) and (|5.10j) . we conclude that for all fixed X £]0,A[, the applications 
z — > Pj(X,X,z) are bounded on K + . Clearly, this result remains true on R by an elementary parity 
argument. Finally, applying the Phragmen-Lindelof's theorem ([7], Thm 1.4.2.) on each quadrant of the 
complex plane, we deduce that z — > Pj{X, A, z) is bounded on C. By Liouville's theorem, we have thus 
obtained : 

Pj(X,X,z) = Pj(X,X,0) , VzeC. (5.11) 

Now, we return to the definition of Pj(X, A, z) for z = 0. We observe first that Fr(x, A, 0) = e tXT x 
and similarly Fr(x, A, 0) = e ,Ar x . This is immediate from the definition of the Jost function. Thus we 



deduce from (|5.5|) that 

P(X,X,0) = e lX ^ x ^ x ^ r \ (5.12) 
Then, putting ([5~l"2l and (IBTTTj) into ([53]) we get 

/ gi(X,X,z) = gi(X,X,z), 

\g 2 (X,X,z) = e - 2 ^ e ^W. 92 (X,A,z), [bA6) 

where we have set 9(X) = X (h(X) - h{X)). 

By Lemma [4~2l the wronskians W(gi 1 g2) — ^(51,52) = i^. Then, a straightforward calculation gives 

e Me(x)-Xc) = l (5 M) 

Thus, by a standard continuity argument, there exists k g Z such that 

6>(X ) = Ac + fcTr , Vl€]0,4 (5.15) 

or equivalently 

^(X) = /i(X) + c+^ , VIe]0,4 (5.16) 
A 



Let us differentiate (|5.16l) with respect to X. We obtain easily 

1 1 



a(h(X)) a(h(X)Y 

and using again (|5.16[) . we have 



(5.17) 



a{x) = a(x + c+ -^) , Viet. (5.18) 
A 

Thus, we have proved the first part of Theorem ll.il 

We are now in position to finish the proof of Thm ll.ll and prove the uniqueness of the mass M, the 
square of the charge Q 2 and the cosmological constant A of the black hole. First, recall that 

Al)= ^)4_^ + ^_|, (6 . 19) 

where r stands for r(x) the inverse of the Regge- Wheeler diffeomorphisn(f|. 



5 We emphasize here that r(x) depends on the parameters we are looking for. 
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To prove the uniqueness of the parameters, we use the following trick. We define the differential 
operator B by : 

g 1 d 2 ^ 20) 

a 2 (x) dx dr' 

dv 

since — = F(r). Using the notation B 2 — B o B, etc., a straightforward calculation gives 
dx 

BW) = + ( , 21) 

9/ 9N 12 M 12 Q 2 , 
5 2 (a 2 ) = 2 — + —£-. (5.22) 

B 3 (a 2 ) = 12M-?i^. (5.23) 

r 

S 4 (a 2 ) = 24 Q 2 . (5.24) 
Now setting x — x + c + and using (|5.18|) , we remark that : 

£ - 1 ^ - 1 d (5 25) 

a 2 (i) c?i a 2 {x) dx 

We apply this differential operator to the equality d{x) 2 = a(x) 2 . To simplify the notation, we set 
f = f(x) and r — r(x). 

Using (|5.24|) and (|5.23|) . we obtain successively : 

Q 2 = Q 2 , (5.26) 
11 M — M . . 

= ^ :=i? - (5 - 27) 

Then, using (|5.22l) . we have 

*£±K E = (M- EQ 2 ) E = ^-M + 2Q 2 E = Q 
2 r 

So, we deduce from (fQ8|) that E = since M, M > 0, i.e. we have obtained M — M and r — f. Using 
now (|5.19p . we get A = A and the proof is complete. 
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